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Abstract 


We study the family of causal double product integrals 



where P and Q are the mutually noncommuting momentum and position Brown¬ 
ian motions of quantum stochastic calculus. The evaluation is motivated heuris- 
tically by approximating the continuous double product by a discrete product in 
which infinitesimals are replaced by finite increments. The latter is in turn ap¬ 
proximated by the second quantisation of a discrete double product of rotation¬ 
like operators in different planes due to a result in HP 15 . The main problem 
solved in this paper is the explicit evaluation of the continuum limit W of the 
latter, and showing that W is a unitary operator. The kernel of W is written 
in terms of Bessel functions, and the evaluation is achieved by working on a 
lattice path model and enumerating linear extensions of related partial order¬ 
ings, where the enumeration turns out to be heavily related to Dyck paths and 
generalisations of Catalan numbers. 

AMS Subject Classification: 81S25, 05A15, 06A07 

Keywords: causal double product, Levy’s stochastic area, position and mo¬ 
mentum Brownian motions, linear extensions, Catalan numbers, Dyck paths. 

1 Introduction 

Following Volterra’s philosophy of product integrals as continuous limits of dis¬ 
crete products [Sla07], quantum stochastic double product integrals of rectan¬ 
gular type have been constructed HP 15 as limits of discrete approximations 
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obtained by replacing stochastic differentials by discrete increments of the cor¬ 
responding processes. Such constructions are partially intuitive in character, in¬ 
volving nonrigorous manipulations of unbounded operators. Nevertheless they 
can be shown to yield explicit rigorously unitary operators which can then be 
shown in some cases HJ12 to satisfy the quantum stochastic differential equa¬ 
tions (qsde’s) whose solutions provide the rigorous definition of the product 
integral. 

In this paper we initiate the much harder problem of constructing so-called 
causal (or triangular) double product integrals in the same way, first construct¬ 
ing discrete approximations by intuitive methods, which are then shown rigor¬ 
ously to enjoy the property of unitarity, which will allow rigorous verification of 
the qsde definitions. 

The Fock space T (7~L) over a Hilbert space T~L is conveniently defined Par92 
as the Hilbert space generated by the exponential vectors e (/) ,/sM, satisfying 

(e (/), e (g)) = exp (/, g) , f,g€U. 

Rectangular product integrals live in the tensor product of two Fock spaces. 
This form of ” double” construction was originally motivated by its use to con¬ 
struct explicit solutions of the quantum Yang-Baxter equation with a given 
classical limit Hud07a , HP05 , of purely algebraic character as formal power 
series. From the analytic point of view, the alternative causal constructs which 
are studied in the present paper which live naturally in a single Fock space are 
of wider interest. 

One example which we consider in some detail is closely related to Levy’s 
stochastic area [Lev51 , and in particular to the Levy area formula for its charac¬ 
teristic function. In effect we replace the planar Brownian motion by a quantum 
version in which the two components are the mutually noncommuting momen¬ 
tum and position Brownian motions P and Q of quantum stochastic calculus 
(CH77 , which despite noncommutativity, can be shown to be independent in a 
certain sense Hudl3 . Other noncommutative analogs of Levy area are based 
on free probability CDM01 ; our own less radically noncommutative form is 
directly related to physical applications HCHJ13 . It may also offer mathe¬ 
matically significant relations, for example to Riemann zeta values through the 
links to Euler and Bernoulli numbers [IT 10 of the classical Levy area formula. 
This is because, while the corresponding probability distribution is the atomic 
one concentrated at zero, it deforms naturally to the classical distribution at 
infinite temperature as the Fock ’’zero temperature” momentum and position 
processes P and Q are deformed through corresponding finite temperature pro¬ 
cesses [CH13 to mutually commuting independent Brownian motions. 

We denote rectangular and causal product integrals by 

II ( 1 + dr )> II ( 1 + dr ) (!) 

[a,b)x[c,d ) <[a,b) 

respectively, where <r a ,b) is the set {(x,y) £ M. 2 : a < x < y < fo} . They are op¬ 
erators in the Hilbert spaces T ( L 2 ([a, b))) ® T (L 2 ([c, d))) and P (L 2 ([a, b ))) 
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respectively. Both are characterised by the generator dr which is a second rank 
tensor over the complex vector space X =C (dP, dQ, dT) of differentials of the 
fundamental stochastic processes P, Q and the time process T of the calcu¬ 
lus. They have rigorous definitions as solutions of either forward or backward 
adapted quantum stochastic differential equations Hudl4 in which b or a in 1 
is the time variable. They are related by the coboundary relation 


II i l + dr ) 

<[a,c) 


(1 + dr) ® I I (1 + dr) 


. <[a,b) 


[a, b) X [6,c) 


I <g> (1 + dr) 

<[b,e) 


in which the Fock space T (L 2 ([a, c))) is canonically split at time b € [ a,c ); 


P(L 2 ([a,c))) 


P(L 2 ([a,b))®L 2 ([b,c))) 

P(L 2 ([a,&)))©P(L 2 ([6,c))) 


allowing it to accommodate the operator ^ ^ (1 + dr). 

A necessary and sufficient condition that they consist of unitary operators 
is Hudl4 that 

dr + dr^ + dr dr' = 0. 

Here the space X =C (dP, dQ, dT) is equipped with the multiplication given by 
the quantum Ito product rule 



dP 

dQ 

dT 

dP 

dT 

-idT 

0 

dQ 

idT 

dT 

0 

dT 

0 

0 

0 


and X © X with the corresponding tensor product multiplication, together with 
the natural involution f derived from the self-adjointness of P, Q and T. 

Two examples of such unitary generators are 

dri = i (dP ® dQ — dQ © dP), 
d7~2 = i (dP ® dP + dQ (g) dQ). 

dri relates to quantum Levy area in which the independent one-dimensional 
component Brownian motions of planar Brownian motion are replaced by P and 
Q. In the same spirit, dr 2 relates to a quantum version of the Bessel process, the 
radial part of planar Brownian motion. The general form of unitary generator 
in which the time differential dT does not appear is [HP15 the real linear 
combination 

dr A , M = -dri + ^dr 2 . 
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In this paper we begin the explicit construction of the unitary causal double 
product integral 


E := JJ (1 + dr\ tll ) 

<[a,b) 

as the second quantisation T (W) of a unitary operator W which differs from 
the identity operator I by an integral operator on the Hilbert space L 2 ([a, b)) 
whose kernel will be found explicitly. 

Acknowledgements: Parts of this work were completed when the Robin Hud¬ 
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Korea, whose warm hospitality is gratefully acknowledged, along with conver¬ 
sations with John Gough, Paul Jones and Janosch Ortmann. Parts of this work 
were completed when Yuchen Pei visited the School of Mathematics of Trinity 
College Dublin, Mathematical Research and Conference Center of the Insti¬ 
tute of Mathematics of the Polish Academy of Sciences, Chungbuk National 
University and Gyeongbokgung. Conversations with Neil O’Connell is also ac¬ 
knowledged. The research of Yuchen Pei is supported by EPSRC grant number 
EP/H023364/1. 

2 The Levy stochastic area 

Before moving on to construct E, let us take a detour and explain the motivation 
of this problem. 

The stochastic Levy area introduced in [Lev51 is defined as the signed area 
formed by connecting the endpoints of a 2-dimensional Brownian path. More 
specifically, it is defined as 

L = \ ( dB h dB 82 - dB 2 si dB l S 2 dsids 2 , 

^ JQ<S 1 <S 2 <t 

where B 1 and B 2 are two independent Brownian motions. The Levy area for¬ 
mula shows the characteristic function of L: 

E n f 1 + y (<*£« dB 2 2 - dB 2 sl dB 2 2 )] = Ee iXL = sechy. (2) 

0<Sl<S 2 <t ' ' 

The Levy area formula has many interesting connotations. For example 
there are connections to integrable systems, Bernoulli and Euler polynomials, 
and hence to the values of the Riemann zeta function [Yor80 . For some re¬ 
cent work and further references see IT10, IT11 . Also, to within normalisation 
and rescaling it is equal to its Fourier transform, the density of the correspond¬ 
ing probability distribution, which is a boundary point of the Meixner family 
fMei34. 

Noncommutative analogues of Levy area have been previously considered 
in free probability CDM01, Ortl3, Vic04 . Also in this connection Deya and 
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Schott [DS13 emphasise the primacy of iterated stochastic integrals which ac¬ 
cords with our philosophy. But in this paper we are concerned with a non- 
commutative analogue of a more conservative kind which arises in quantum 
stochastic calculus [HP84, Par92 , regarded as a noncommutative extension, 
rather than a radically noncommutative analogue, of Ito calculus. This allows 
a very natural variant of the area to be constructed using the minimal one¬ 
dimensional version of the calculus. It may be regarded as a response to the 
call ApplO for a study in this quantum context of some of the deeper prop¬ 
erties of Brownian motion, as well as a furtherance of the theory of quantum 
stochastic product integrals [Hud07a ( Hud07b, HP81 . 

By replacing B 1 and B 2 with P and Q, the iterated quantum stochastic 
integral 


K(t) = 7 : [ ( dP x dQ y - dQ x dP v ) 

has some interesting properties [Hudl3, CH13 . For example it is evidently 
invariant under gauge transformations, which replace ( P,Q ) by (P 8 ,Q 8 ) where 

P e = P cos 0 — Q sin 9, Q e = P sin 0 + Q cos 9; 

equivalently the corresponding creation and annihilation processes are multi¬ 
plied by e ±l8 . In particular, taking 9 = — ^ it is invariant under the replace¬ 
ment ( P,Q ) by ( Q,—P ). Thus, unlike the separate processes P and Q , it can 
be canonically “rolled” onto a (one-dimensional) Riemannian manifold, and its 
multidimensional version FV10 can similarly be rolled onto a multidimensional 
manifold, with possible applications to quantum stochastic proofs of index the¬ 
orems, by identifying the canonical Brownian motion on the manifold generated 
by the Laplacian as P 8 with arbitrarily chosen 9. 

It can also be verified CH13 that all moments of K(t) vanishes in the 
vacuum state, so that K(t ) vanishes in a probabilistic sense, even though it is 
not the zero operator. 

But it is not K which is the main object of study. Because exp(a + b) ^ 
exp a exp & when a and b do not commute, the exponential 

exp(i\K(t)) = exp ( ^ f ( dP x dQ y - dQ x dP y )\ 

does not reflect in a coherent way the continuous tensor product structure under¬ 
lying the quantum stochastic calculus. Thus, motivated by the hope of finding 
quantum extensions of, in particular, the Levy area formula (2), and associated 
relations with Euler and Bernoulli polynomials (IT11 we investigate the double 
product integral 


JJ(l + dri)= (l +y( dPxdQy - dQ x dP y )\ 

<[a,b) a<x<y<b v ' 
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However, as it turns out, the more general object 

E = JJ (1 + drx tfl ) 

<[a,b) 

= f 1 + —{dP X dQy — dQxdPy) + —(dPxdPy + dQxdQy ) J . 

a<x<y<b ^ ' 

is more fundamental and, surprisingly, simpler to study. 

3 A discrete double product of unitary matrices 

The first stage of the construction of E is similar to that of the rectangular case 
construction outlined in [HP15 , in that we approximate II (1 + dr\ by 

a discrete double product ni<j<fc<Ar (j + r\,^j , where S^ k r\^ is obtained 

from dr\ by replacing each basic differential dX £ { dP , dQ} contributing to 
dr £ X g) X in the first copy of X by the j- th increment X Xj — Xj.and in 
the second copy of X by the fc-th increment X Xk — X Xk _ 1 over the equipartition 



Thus, for example, 


r l — 2 ((Aj Pxj-i) < 8 > ( Qx k Qx k _ i) 

— (Qxj — Qxj- 1) ® {Px k — Px k _ 1 )) ■ 


Introducing the standard canonical pairs (pj,qj) ,j = 1,2,..., N, given by 



which satisfy the canonical commutation relations 

[Pj,Qk] = ~2i5j,k, [ Pj,Pk ] = [qj,Qk\ = 0, 


(3) 


we write 



and more generally 



Our approximation is thus 



( 4 ) 
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for large N. 

Temporarily let us fix j < k and write (p, q) = (pj,qj ), {p' ,q') = ( Pk,Qk ) so 
that 

\p,Q] = -2 i, [p',q'} = -2 i, [p, q'] = [q,p} = \p,p') = [q,q'] = 0. (5) 

We recall [Par92 that, for an arbitrary Hilbert space 7~L and vector f £ 7~L the 
corresponding Weyl operator W (/) is the unique unitary operator on T (7~L) 
which acts on each exponential vector e (g), g £ H as 

W {f) e (g) = e -5ll/ll 2 -U,g) e (f + g). 

The Weyl operators satisfy the Weyl relation 

W(f)W(g) = W (/ + g). (6) 

A convenient rigorous realisation of two canonical pairs satisfying the commu¬ 
tation relations (5) can be constructed in terms of the one-parameter unitary 
groups of which they are the self-adjoint infinitesimal generators, which are 
Weyl operators on the Fock space T (C 2 ) over C 2 . Regarding C 2 as a space of 
column vectors, we take 

e ixp = W (Or, 0) r ), e ixq = W ((-ix, 0) T ), 
e ixp ' = W({Q,x) T ),e ixq ' = W ((0,-*x) r ). 

for arbitrary x £ R, noting that these four families of Weyl operators are indeed 
one-parameter unitary groups, and that the commutation relations (5) follow 
by parametric differentiation, for example from the relations 

W((z,0) T )W(Hp,0) r ) = e 2ixy W ((—iy, 0) T ) W ((x, 0) T ), 
W((0,x) T )W((0 : -iy) T ) = e Mx yW((0,-iy) T )W((0,x) T ), 

all of which are consequences of (6). 

Theorem 1 below, which is proved in [HP15 , gives a corresponding rigorous 
explicit form of the self-adjoint operator 

L (A, p) = A (pq' - qp') + p (pp' + qq') 

in this realisation. Before stating it we recall [Par92 that the second quantisa¬ 
tion of a unitary operator U on a Hilbert space T-L is the unique unitary operator 
r (17) on T (TL) which acts on the exponential vectors as 

P (U) e(f) = e (Uf). 

It is related to the Weyl operators by 

T(U)W(f) = W(Uf)T(U) (7) 

for arbitrary f £71. Second quantisation is multiplicative, in the sense that 

T(U 1 U 2 )=T(U 1 )T(U 2 ) (8) 

for arbitrary unitary U\,U 2 . 
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Theorem 1. L( A,/r) generates the one-parameter unitary group 


e ixL(\,ii) __ p 


cos (2x |^|) —e sin (2x\u\) 

e 1 ^ sin (2x \v\) cos(2a’|^|) 


x £ R 


where v = A + ip, = e 1 ^ |^| and the matrix operates on column vectors in C 2 by 
multiplication on the left. 

We now use Theorem 1 to construct an explicit second quantisation of the 
approximation (4). 

Let us first construct a different realisation of the canonical pairs (pj, qj ), j = 
1, 2, n, satisfying (3) in the Fock space T (C ra ) over C n , by defining 


= W {x£j) , e ixq > = W {-ixej) 


where (e^)” =1 is the standard orthonormal basis of C n , £j = |^0,1 ,0,0 
Correspondingly, in view of Theorem 1, each operator 


ex P ( A “ 9jPk) + P (PjPk + qjQk))j 

is realised as the second quantisation T(R^ k ) where 


jdN _ 

n j t k ■ — 


U) o 


(k) o 


U) 

o 


COS 




H 


s in(^^|i/|) 


(k) 

0 


H sin (VM) 


COS 


(Vl-l) 


0 

0 

0 

1 


In view of the multiplicativity property (8) the discrete double product (4 ) is 
correspondingly realised as the second quantisation of the product 

n R xk- (9) 

1 <j<k<N 

We now embed the matrix (9) as a unitary operator Wn on L 2 ([a, &[) by map¬ 
ping the standard basis of C N to the orthonormal family (xi, X 2 ; ■ Xjv) of 
normalized indicator functions 



[*i- 


Xi ( X ) = 

















By definition Wn acts as the identity operator I on (yi, X 2 , ■■■, Xn) ± ■ 

Our objective in the remainder of this paper is to find an explicit form for 
the (weak) limit 

W = lim W N 

N—t OO 

and to prove that R is unitarity. The corresponding problems for rectangu¬ 
lar unitary product integrals was solved in outline in (HP15 . The causal case 
considered here is considerably more difficult, because the method of iterated 
limits which reduces the rectangular case to a double application of the time- 
orthogonal unitary dilation of HIP82 , is not applicable. Instead a combinato¬ 
rial argument based on a lattice path model is used. For a similar alternative 
approach, avoiding the iterated limit technique, to the rectangular product in 
the particular case of the generator dri corresponding to the quantum Levy 
area, see [HJ12J; however the combinatorics for the rectangular case is much 
simpler than here and it has no direct relation to Levy area. 


4 A lattice path model and linear extensions of 
partial orderings 

So we want to calculate the limit of the triangular double product of IV x N 
matrices 

n Kk- (io) 

1 <j<k<N 

Here, for elements xj t k of an associative algebra having the property that Xj } k 
commutes with Xj>k' whenever both j ^ j' and k ^ k! we define the ordered 
double product rii<j<fc<Ar x iM by any of the equivalent prescriptions 


n 

1 <J<k<N 


N—l 

x j,k = 

i=i 


N 

N 

k -1 

n x x k 

=n 

n x j,k 

_k=j +1 

k =2 

i=! 


ijV(JV-l) 

n x jr,k r 


r—1 


where ^(ji, fei) , {j 2 ,k 2 ) , (jijv(jv-i)) fc iJV(JV-i))) is any ordering of the \N(N- 
1) pairs (j, k), 1 < j < k < N which is allowed, meaning that it has the property 


that 

(j r , k r ) precedes (j s , k s ) if both j r < j s and k r < k s . (11) 

In constructing the limit as N —> oo we use the small angle approximations 
for sine and cosine, so that 


C0S TrM -fsin^M 
fsin^l v\ cos^H 


= 1 


b — a (0 — 


N 


v 0 


0{N~ Z ) 


hence 


R xk = 1 + ^ K T a \ - v I Xj) (Xk\ + v | Xk) (Xil) + 0(N 2 ). 


N 
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When there is no ambiguity, for any integers j and k, we use abbreviations 
| j) := | Xj) and (fc| := (xfe|- Then the product (10) becomes 


Wn ^ 


n u+ 

1 <j<k<N V 


jb-a) 

N 



=: W N 


where 

Z N {j, k) = -v | j) (k\ + v |A:) {j\ . 

To compute this, we introduce and work on a lattice path model. Consider 
a lattice L s := {(to, n) : 1 < m < s, 0 < n < 1}. We call (to, l)i< m <s the upper 
vertices, and (TO,0)i< m < s the lower vertices. Denote by II S the set of lattice 
path n = (rrii,bi)l =1 satisfying the following two conditions: 

1 . to* = * for i = 1 ,..., s 

2. there does not exist an i such that bi = 6,;_|_i = 0 

For convenience, we write n(i) = bi and let 7r = (n{i))i- We call any tt £ II S a 
path of length s — 1. It is straightforward to verify by induction that 


I n s | — Fib s +2 — 


$ s + 2 — (—<f>) —s—2 

7 ! 


where Fib ra is the nth Fibonacci number and $ is the golden ratio v/ ^ +1 . 

If we assign weight 9(v) to each vertex v in L Sl then we can define the weight 
0(n) of a path n £ II S by the product of the weights of its vertices: 


6 >( tt ) := Y[d(i,Tr(i)). 

j=i 


For any s-array of pairs {pij : 1 < * < s, 1 < j < 2}, define its associated 
weight 6 p (v ) for any v = (to, b) £ A to be 


0 p {v) 


v\pi 2 ){pn\, if b = 0; 

—v \pn) {Pi 2 \ , if 6=1. 


Finally, define the weight 0 p (ir) of a path in the same way as before. 

For example, if we label the vertices by their weights associated to p , then 
the following is a path of II5: 


— v \vil) (P12I- v |P2l) (P22I — v |p3i) (P32I — v |p4,l) (P4,2 1- v |p5,l) (P5,2 



+ V |i>12> (pil | +" IP22) (P21 1 + V |j> 32 > (P 311 +" 1^4,2) (p 4 ,l | +V |P 5 , 2 ) (P 5 ,l 

but not the following because the third edge connects two bottom vertices: 
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-HPn) (P12I 


|P2l) (P22| -V |p 3 i) (p 32 | |P4,l) (P4,21 |P5,l) (P5,2 


+ HP 12 ) (Pll| +^|P22)(P2l| +^|P32) (p 3 l|-MP4, 2 } <JP4,l| +^1^5,2) (P5,l 

Any s-array of pairs p = (Pij)i<i< s ,i<j <2 satisfying the following condition 


-^Pi,i=Pi+i,i d" Ip», 2 =p«+ 1 , 2 d~ ^pi,2=pi+i,i 


Pi, 1 7 ^ Pi+ 1,2 


( 12 ) 


can be associated with a path tt p £ II S in the following way: 

{ (0,1), if Pi,i = Pi+ 1,1 
(1,1), if Pi ,2 = Pi+1,1 ■ 

(1,0), if Pi ,2 = Pi+ 1,2 

Note that this is equivalent to 

S 

0p(7fp) = II ZN (PiA’Pi,2)- 


Lemma 2. 


where 


i=i 


N(N—l)/2 

W N = I + ^ w 3> jy, 

S—1 


w s,n = 1 ) y^flpM) 


IV 


(*) 


where the domain (*) of the summation is 

(1 < Pi,j < N,YL < i < s, 1 < j < 2) AND 
{jpi, 1 < Pi+1,1 < Pi ,2 = Pi+1,2 0-R 
Pi,i < Pi, 2 = Pi+ 1,1 < Pi+1,2 OR 
Pi, 1 = Pi+1,1 < Pi ,2 < Pi+ 1 , 2 , VI < i < s - 1) 

Proof. For any rearrangement (.7,;, fci) 1 < . < w(jv-i) of {1 < j < k < N} satisfying 

fin. 

N(N — l)/2 N(N-l)/2 s 

n (i+z(j,k))= n (i + z( Jl ,k i ))=i+ Y, En Z(j?rl , Prf ), 

1 <j<k<N i= 1 s=l (**) r=l 

where domain (**) is 


(Pll, Pl2 ) , (P21, P 22 ), ■ • • , (Psl,Ps2) 

is a subsequence of (ji, fci), (j 2 , k 2 ),(j N (N- 1 )/ 2 , fcjv(jv-i)/ 2 )- 
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Now for the product J~[r=i Z{p r i,p r 2) to be nonzero, the product of P- 1 , 2 ) x 

Z(pi + i t i,pi + i^) for each i has to be nonzero, that is 

(- \Pi, 1) (PiM + \Pi, 2) <P»,i|)(- |Pi+i,i> {Pi+iM + \Pi+ 1,2) {Pi+ i,i|) ^ 0 . 

This in turn is equivalent to 


(Pi ,2 = Pi+1,2) V (p ij2 = Pi+1,1) V (pi, 1 = Pi+1,1) V (pi,i = Pi+i, 2 ). 


We analyse these four possibilities one by one. 

1 - If Pi, 2 = Pi+1,2, then by ( 11 ), and since (pi,i,Pi, 2 ) 7^ (Pi+i,i,Pi+i,2), only 
whenp^i < Pi+ 2 ,i can the product be nonzero. In this case the coordinates 
are ordered as p t + < Pi+1,1 < Pi, 2 = Pi+i,2- 

2. If Pi, 2 = Pi+1,1, then since p f) 1 < p» )2 = Pi+1,1 and p» )2 = Pi+1,1 < Pi+1,2, 
we have that (11) is satisfied. Therefore this case is also included / per¬ 
mitted in the product. The ordering of the coordinates is pi+ < p^ j2 = 
Pi+ 1,1 < Pi+ 1 , 2 - 

3. If pi+ = Pi+ 1 , 1 , then similar to Case 1, the coordinates have to satisfy 
Pi, 1 = Pi+1,1 < Pi, 2 < Pi+1,2 for the product to be nonzero. 

4. If p^ 1 = Pi+i, 2 , then p i: 1 = p i+ i, 2 > p +1 and p ii2 > p+i = Pi+1,2 violates 
(11), hence this case never happens. 


The three feasible cases are illustrated as below. 

i i + 1 i i + 1 i i + 1 

/ 

Pil <Pi+l,l Pil <Pi+l,l Pil = Pi+ 1,1 

A A A A- A A A 

Pi 2 =Pi+l ,2 Pi 2 <Pi+l ,2 Pi 2 <Pi+l ,2 



(13) 


The concatenation of these edges gives a path in II S . Case 4 corresponds to 
a horizontal bottom edge in the path which is not allowed in the definition of 
U s . Therefore we have established a correspondence between the possibilities of 
orderings in the product and II S . □ 


Denote by ^4* the set of s-array pairs p satisfying condition (*), and f:= 
{p € A* : tt p = 7 r}. Then 


W S ,N 



Given a path n £ II S , by the correspondence in (13) there exist mi, m 2 , ..., 
m s -1 S {1, 2} such that for anyp G fi*-, p x ,m x = P x +i,m' , where rn' x := 3 —m x . 
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Therefore, is characterised by a partial ordering on the s + 1 coordinates 



nates of p. This also shows we can associate tt with (mi,TO 2 ,... ,to s _i). In 
the following we do not differentiate between tt and the corresponding partial 
ordering. 

Any ordering n £ II S can be decomposed into (strict) total orderings without 
any repetition of the essential coordinates and those with repeated essential 
coordinates. We call any the former orderings B a linear extension of 7 r which is 
denoted by B b tt, and the latter degenerate orderings, which, for reasons that 
will emerge in the proof of Lemma 3 are ignored. Thus we have 


fin- = |_| A? U set of degenerate orderings. 

Bhn 


For any p £ fl n , there exists a B b n such that p € B. Denote (jp,k p ) = 

(j)- 1 11 the total ordering imposed by B, let tb be the number of 
essential coordinates less than j p and r' B the number of those greater than k p . 
That is, the essential coordinates are ordered as follows, 

1 < h < ■ ■ ■ < l rB < j P < nii < ■ ■ ■ < m s _i_ rB _ r ' B < k p < ni < ■ ■ ■ < n r ' B < TV, 


if r B +r' B < s 


1 < /l < * * * l s — r ' B <C k p <C TTl\ ‘s * * * 's TTl rB J rr ' B — s —\ < j p Tl 1 * * * < n s — rB ^ A^, 


if r B +r' B > s 


We call ( r B 'b ) the rank of B. 

Let e( 7 r) be the number of upper vertices of the path tt. Since horizontal 
edges between lower vertices are not allowed, there is at least one upper vertex 
between two consecutive lower vertices, hence 


2 e( 7 r) > s — 1. 


The location (upper or lower) of the first vertex of n, the number of upper 
vertices e( 7 r) and the parity of the length of n together determine the number 
of horizontal edges in tt. The cases when e(- 7 r) ft: are “saturated”, meaning 

there is no horizontal edge n. This will be later specified and exploited in the 
proof of Lemma 12. 

The weight of tt is 


0 p(n) = (-9yWv s -^ | j p ) (k p \. 


So 
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where f s ,jv is the contribution from the degenerate orderings, on which one can 
carry out the same calculation for w s _m below, and that 




\ N Ei <i l<---<lr<j<'mi<---<m s _ 1 _ r _ r / <fe<ni<-<n r / <N I i> (fc| , 

(~aT i ) Ei<; 1 <-"<l s _ r / <k<m\ <---<m r _|_ r /_ s _ 1 <j<ni <---<n s _ r <N I j) ( k \ , 


r + r' < s 
r + r' > s 


For example, for the following path 7 r of length 2, 



Pl,l = P2,1 <P 3,1 
AAA 
Pi,2 < P2,2 = P 3,2 


The ordering of the essential coordinates imposed by 7r is: 

(pi,l < Pi,2 < P2,2) A (p 2 ,l < P 3 ,l < P2,2), 

and the non-repeated starting and ending coordinates are j p = p±2 and k p = P31 . 

The total ordering decomposition of is 

= {pil < P12 < P 31 < P22} LI {pn < P31 < P12 < P22} U {pn < P12 = P 31 < P22}- 


The last term is a degenerate case as p \2 is repeated. There is only one upper 
vertex, hence this path contributes — v\v\ 2 {Hff (2, 3) + H^( 3,2)) to in s ,jv- 
Define the Volterra-type kernels > b a ( x,y ) := l a < y <x<b and < b a (x,y) :=> b 
(y,x), and [■ m,n,p\(x,y ) := {v ~ x) (b ~f } and [m, n,p] f (a:, y) := [m,n,p](y,x). 

The asymptotics of if™ can be written down explicitly. 


Lemma 3. Hj?(r,r') converges weakly to an integral operator H s (r,r'), with 
the integral kernel 


hs (r, r') 


[r,s —1 — r-r',r'] < b , r + r'<s 

[s — r',r + r' — s—, s — r] f > b a , r + r'> s 


Proof. Suppose r + r' < s (the case r + r' > s can be done in the same way). 
Then 


Hg (r,r) = 


b — a 
N 

b — a 
N 


H H 1 

l<j<k<N l<l 1 <---<l r <j<mi<---<m s _ 1 _ r / _ r <k<ni<---<n r /<N 

z ' \ r J \s — 1 — r — r t 

1 <j<k<N \ / \ / 
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We denote Ajv := ^2, then the kernel of H^(r,r') is 


hs (a ’-,y) 


E 1a i ( x ) Ia * (y) 

1 <j<k<N 


1 

r\(s - l - r - r')\r'\ 


r —1 s—2—1 — r' r'— 1 

x (Xj-i - a - aAjv) (^ fc -i - Xj - /3A N ) (b - x k - tA n ). 

ck=0 /3—0 7=0 


This, as A —> oo, converges weakly (as an integral kernel) to [r, s — 1 — r — 

r',r'] < b a (x,y). □ 

It can also be seen from the proof of this lemma that the degenerate order¬ 
ings contribute 0 to the total sum. More specifically, the degenerate version of 
h^ ( x , y) where there are d repeated essential coordinates is 




r— 1 


E 

1 <j<k<N 


1 a, {x)i Ak (v)- 


!(s - 1 - r - r')b 


s—d— 2 —r—r' 


(i b-a) d 
N d 


r' — 1 


x ]](i H -a-aAiv) {x k -i - Xj - /3A N ) (6- x k - yA^) -t 0 

a=0 /3=0 7=0 


as N —> oo. We will examine carefully the rate of convergence of this lemma and 
the (in)significance of the degenerate orderings later in the proof of Theorem 7. 
This lemma immediately gives the following corollary. 

Corollary 4. There exist two integer arrays ( Dm,n,p-, q )m,n,p>o,o<q<m+n+p+i 
and (Em,n,p;q)m,n,p>o,o<q<m+n+p+i such that w Sj n converges weakly as N — > oo 
to w s with kernel 


fs <a +ffs >a 


where f s and g s are defined by 

S 

fs= E '52 D ™„n,p-,q(-v) q v s ~ q [m,n,p\, 

m,n } p>0,m-\-n-\-p=s— 1 q—0 
s 

9s= E E E m,n,p-,q{-v) q V S ~ q [m, n,p]f 

/ m,n,p>0,m-\-n-\-p=s—l q =0 

Indeed, D mtn ^ p . q (resp. E m ^ n ^ q ) enumerates the linear extensions of all 
possible paths of length m + n+p with q upper vertices and rank ( m,p ) (resp. 
(■m + n+ 1 ,n + p+ 1)). 

Corollary 5. The functions f s and g s both are symmetric in the following 
sense: 

fs{x, y) = f s {a + b — y,a + b — x), g s (x,y) = g s {a + b - y, a + b - x) 
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Proof. This follows from the fact that the path inversion (*, bi ) H > ( i , is 

a weight-preserving bijection between II S and itself. □ 


For example, some calculation yields 

fi = -^[ 0 , 0 , 0 ], 

h = -k| 2 [0, o, 1] - M 2 [1,0,0] + v 2 [0,1,0], 
h = (v\v\ 2 - ^ 3 )[0, 2, 0] + v\v\ 2 [0,l, 1 ] + I'M 2 [1,1,0] - ^|^| 2 [1,0,1], 

91 = -^[ 0 , 0 , 0 ] f , 

92 = 0 , 

93 = -i/|z/| 2 [l,0,l] t . 

The following three theorems are the main results of this paper: 

Theorem 6 . The closed form expression of D and E are: 

f («-i) — (")» 2 q>m + n + p 

D m ,n, P -,q = l ( q " m ) - Q), 2q = m + n + p (15) 

[0, 2q < m + n + p 

Em,n : p;q — ^-m=p=q,n =0 


Proof. See Section 5, □ 

Theorem 7. The operator Wn converges weakly to 

W = 1 + ^2 w s 

s> 1 


Proof. See Section 6 , □ 

For j > 0, let Bj be power series in two variables related to the Bessel 
functions of the first kind Jj. 

Bj(x,y) := Y, ( l ( n+ X j)\ n \ V = 

n> 0 ' 

Let I be the identity, then the kernel of the operator W — I can be written in 
terms of Bj. 
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Theorem 8. The integral operator W — I has kernel 


ker(W - /)( x, y) = [vB 0 ({y - a)\u\, (b - x)\v\) + \u\B 1 ((b - a)\v\, {y - x)\v\) 

- (" + - x )W\A b -a)W\) (j^j j < b a(x,y) 

+ uB 0 ((y - a)\v\, (b - x)\u\) > b a (x,y). 

Moreover, W is unitary. 

Proof. See section 7 □ 

For example, when /i = 0 and A > 0, the kernel of the operator corresponding 
to the Levy stochastic area is 


ker(TT — I)(x, y) 


^AB 0 ((j/ - a)A, (b - a;)A) + XBi((b - a) A, (y - x)\) 


- 2 A 

q> 0 


x)\,(b 



< b a ( x ,y) + ^ B o{{y 


a) A ,(b 


x )*) >a (x,y)- 


Moreover by plugging D mt n tP . q and E m , n ^ q into the integral identity (18) below, 
the unitarity of W implies the following combinatorial identity: 


Lla,/3,7;£ ^a=7=£—1,/3=0 


P+ 7 +l=a=C 


'a + 7 — 

7 

ol 7— OL-\-m at-\-/ 3 — 7—ra+p— 1 

^ ^ ^ ^ ^ ^ -^ra,n,a+/ 3 —7—m—n+2p—1;£—7+p— 1 

m =0 p —0 n =0 

' ^7 — a + ra + n — p\ f 7 ^ 

n y Vp 

a- /3 7—1 7—1—ni 7— 1—n-i—ri2 £ 

EEE E E E(-i ) a+7_roi -" l - pi+ro2 

mi—0 7712 —O ni=0 n. 2 —0 pi=0 ti=0 

-^mi,ni+/ 3 —rri2,pi ;ii ^m2+a-mi ,ri2,7—1—ni —n,2—pi ;a+7—mi—ni —pi+m.2 +ti 


a 

raw \ra 2 


ni + n 2 \ 7 — 1 — ni — n 2 


ni 


Pi 


= 0 . 


5 Dyck paths and Catalan numbers 


For to e Z>o and n£Z, dehne the binomial coefficient the usual way 



to ! 

77 .! (to. — 
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For integers m,n,p define a double generalisation of the Catalan numbers and 
the Catalan’s triangle 


C 



m + n 


m 


( to + n 
to + p + 1 


For a,m,n,p € Z>o, denote by T aim>njP the set of lattice paths {pijffl o" suc ^ 
that p 0 = a, |Pi - Pi-i | = 1, Pi > -p, Pm+n = a + m - n. That is, T a ,m,n, P 
is the set of Dyck paths starting from a, having to up-steps, n down-steps that 
never cross the line y = —p. By the reflection principle we obtain the following 
lemma, which shows these numbers have a similar combinatorial interpretation 
to the Catalan numbers. 

Lemma 9. When m,n,p> 0 and m — n> —p — 1, C m ^ p = |To )mi „ jP | . 

The doubly generalised Catalan numbers have been discussed in e.g. (Reul4 . 
When to, n > 0 and p = 0, C mi „ : o is reduced to the (to, n)th entry in the 
Catalan triangle (OEIS:A009766) which we denote by C m y, furthermore when 
to = n, C ntn is the nth Catalan number which we denote by C n . 

The following recurrence relation will be useful: 

Lemma 10. If n> 0, m > p and m + n+p+l>0, then 



Proof. We first show a basic version of this formula is true: for n > 0, to > p > 0, 


p 



This can be proved using a combinatorial argument similar to one used to prove 
the recurrence relation of the Catalan numbers which is a special case of the 
identity above: 


p 



Define a “stopping time” a on To irra +„+i iP; o by 

cr(p) = max{* > 0 : pi = n}. 
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Then 


P 

Cm+n+l,p = ^ ^ \{p ^ ^0,m+n+l,p,0 • &(p) = H~ 7T-}| 
k—0 

V 

~ \ T 0,k+n,k,0\\T n +l ,m-k,p—k,n -\-11 

/c=0 

P V 

= ^|T 0 , fc+ „ iM ||T 0 ,m—fe,p— fe,01 — ^ ^ c m - k , p - k . 

k= 0 fc=0 


If the condition n > 0, m > p are retained, but p < 0 and m+n+l+p> 0, 
then the LHS is zero because the domain of the summation is empty. The RHS 
is also zero because = (”«««) = 0. 

Since m > p, we have p < L^r^J ■ Moreover, for any k £ ( p , 
Cm-k, P -k = ( m ^fc 2fc ) — (™-V+T) = Therefore we can extend the domain of 
the summation from 0 < k < p to 0 < k < . □ 

Lemma 11. For any B b 7 r £ II S; if n(0) = 1 f/ien = 0, and if tt(0) = 0 
then r\B > 0. If 7r(s) = 1 f/ien = 0, and if n(s) = 0 then r' B > 0. 

Proof. We show the claim for 7 r( 0 ), as the one for 7 r(s) can be deduced from 
the symmetry property. If 7 r( 0 ) = 1, then 7 r(l) = 0 or 1. If 7 r(l) = 0 then 
by the correspondence (13), for any p £ i the first four coordinates have 
the ordering j p = pip < p 2 ,i < Pi ,2 = P 2 , 2 - Since pip < Pk,i,k > 2, and 
pi,i < Pi ,2 < Pk, 2 ,k > 2. Thus j p is the smallest (essential) coordinate and 
rs = 0. If 7 r(l) = 1 then j p = pi,i < pp 2 = P 2 ,i < P 2,2 hence it’s also the 
smallest coordinates and re = 0 . 

If 7r(l) = 0, then 7r(l) = 1 and by (13), for any p £ fl n , the first four 
coordinates are ordered as p ip = p 2j i < pi, 2 = j P < P 2 , 2 - Hence j p is greater 
than at least one other essential coordinate and r B > 0. □ 

In some extreme cases the coefficient ? can be calculated directly. We 

denote by Dm,np;q the contribution to D mt n tP . q from path n. 

k 

Lemma 12. • (Case A) Do^ 2k ,o-,k+iDo 2 fe,o-fc+i = Ok- Conversely, if m = 

0 , 2 q = n + p + 2 , then Dm^p-q > 0 only if p = 0 , 2 q — 2 = n. 

• (CdSC Bj DQ n 2k— n+l;fc+l — ^ 0 ,n, 2 fc—n+l,fc+l — ^k^n—k for 0 ^ Ti ^ 2k~\~ 1. 

A fc / 

• (CdSC C) D2k— n+l,n, 0 ;fc+l ^2k-n-\-l,n,k-\-l — ^k,n—k for 0 ^ 71 ^ 2 k ~\~ 1. 

A k 

• (CttSe D) D r ,2k—i — r',r';k = -^r,2k—r—r',r'’k = ^k—r,k—r',r —1 = Ck—r',k—r,r' — 1 
for 0 < r + r' < 2 k. 

Moreover, E rnn p-q = I rn _p_g ?n _Q. 
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Proof. First we show the first identity in each case. In Case D, for there to be 
fc upper vertices and k + 1 lower vertices, the path can only be A k . 

For Case A, since the rank is (0,0), by Lemma 11, any path 7 r contributing 
to I?o,2fc,0;fc+i has to begin and end with upper vertices. Removing these two 
vertices resulting a path of length 2 k — 2 , k — 1 upper vertices and k lower 
vertices, which is the same as Case D. Therefore n = V fc . 

With the same arguments the paths for Case B and C are also determined 
to be \A k and A k / respectively. 

Now we show the second identity in Case D, as Cases A, B and C are simpler 
variations of D and can be verified similarly. We achieve this by associating par¬ 
tial orderings with sets of Dyck paths. The path n = A k imposes the following 
ordering of the essential coordinates: 


Pl,l 

A 

< 

P3,l 

A 

< 

P 5 ,l 

A 

< . . 

. < 

P 2 k-l,l 

A 

< 

k p — P2fc+i,i 
A 

j P = Pi, 2 

< 

P 2,2 

< 

PA ,2 

< . . 

. < 

P 2 k- 2,2 

< 

P 2 k ,2 


We relabel these coordinates by tip := p+i, £ 2,1 = P 3 ,i and so on, to obtain 

1 1.1 < £2,1 < £3,1 < ■ • ■ < tk, 1 < ffc+ 1,1 

A A A ... A A (16) 

1 1.2 < £2,2 < £3,2 < < tk , 2 < tfc+1,2 

There is a one-one correspondence between the set of all linear extensions of 
this partial ordering (namely {B : B b 7r}) and To,fc+i,fc+i,o- The Dyck path 
p corresponding to the linear extension < t m2 ,b 2 < ••■ < t mk+1 ,b k+1 is 

defined by 

/^(* _1 ) + 1 > if = 1 
' ,(,) = \,(i-i)-i, if h = 2 

Clearly, the rank of a linear extension B being (r,r') is equivalent to the cor¬ 
responding Dyck path starting with r up-steps followed by a down-step and 
concluding with one down-step with r' up-steps. These cut off the first r + 1 
and the last r' + 1 steps corresponds to T r -i,k-r,k-r’, o- Therefore 

A fc 

D r , 2 k-i —r',r';fc = |T) — l,fc-r,fc-r',o| = \Tq, k-r,k-r' ,r-l\ = &k-r,k-r' ,1—1- 

If r = 0, r' = 0, then by Lemma 11 D^ 2k _ r _ r , r ,. q = 0, which agrees with 
Ck- r ,k-r',r- 1 as well. On the other hand, since the paths of Xo.fc.fc.o only have 
k up- and down-steps, the LHS is 0 if r > fc or r' > fc, which agrees with the 
right hand side. 

Finally, the A k are the only possible paths to contribute to the coefficients E, 
which record the instances when k p < j p . The corresponding linear extension 
B h A k is < f 2 ,i < ••• < ffc+ 1,1 < fi ,2 < * 2,2 < ••• < tfc+ 1 , 2 }• For any 
other paths, by Lemma 11, any path starting with \ or ending with / has j p as 
the smallest or k p as the greatest essential coordinate; on the other hand, any 
horizontal edge will result in j p < k p . □ 
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From the above proof we can deduce a stronger version of Lemma 11: 
-Dq np . q / 0 only if 7r(0) = 1, and for to > 1, -Dm,n, P ; 9 7^ 0 only if 7 r( 0 ) = 0 and 
begins with A" 1 ’ 1 . We also refer to this stronger version as Lemma 11 

Proof of Theorem 6. Case D covers the 2 q = m + n + p case; moreover, by the 
same argument as in the proof of the first identities in each case of Lemma 12, 
for 7 r £ I l !n , 


e(7r) + l>TO + n + p+ l — e(7r) 

therefore q 7 ^ 0 only if 2q > m + n + p. Thus it suffices to show 

Dm,n,p;q — ^P TO T 77 . T p. 

We group the paths into ones starting with A k /— (call the set of such paths 
n Afc ) and ones starting with V fe — (call the set of such paths II vfc ). Then by 
Lemma 11 -D m ,n, P ;< ? are contributed from II Afc if k — 1 > to > 0, and from II V 
if to = 0 : 


D 


m,n,p;q 



k>m— 1 X/7rsn Afc ^m,n,p;g5 
k >0 ^7rCll vfc 


TO > 0 
TO = 0 


where n is the contribution from path 7 r to the coefficient D m ^ n ^ q . There¬ 
fore we divide the proof into two cases, to = 0 and to > 0. The formula for 
D m ,n,p;q with to + n + p < 2 can be verified by hand (the reader can check their 
calculation against (14)), so we assume the formula is true for s < m + n + p, 
and we want to use induction to verify the formula of -D m ,n,p;g in general. 


5.1 m > 0 

When to > 0, for any 7 r £ II Afc , 

7 )7T _ 7-. 7-\#2fc + 2TT 

““ i 'm,2fc+l-m,0;fc-|-l- LA o,m-|-ra-2fc-2,p;<j-fc-l > 

where 9 r : II S —> II s _ r Vs > r is the shifting operator such that (■ 0 r n)(j ) = 
t r(j + r). 
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Summing over k and it £ II Afc we have 

Dm,n,p\q = ^ 1 ^ 1 D m ,n,p-,q 

k>m— 1 7ren Afc 


n n @ 2 fe+ 2 ir 

/ , iy m,2fc+l-m,0;fc+l-^o,r7i+n-2fc-2,p;<j-fc-l 

k=m— 1 7rGlI Afe 

[ n 1± g = 2 i 

= ^ ] ^m,2k+l — m,0- : k+1^0,m+n—2k—2 1 p\q—k—l 

k—m—l 
L m4j^2j 

y \ ^k,k-m+l^-'q—k—2^m+n—q—k 

k=m—l 

L n-mj 

— y \ Ck-\-m—l,kCq—m—l—k,n—q-\-l — k 
fc=0 


where the last two equalities comes from the Case C and the induction assump¬ 
tion. 

To apply Lemma 10, we check the three conditions hold: (1) The condition 
“n > 0 ” becomes m — 1 > 0 : this is correct as m > 0 . ( 2 ) “to > p" is 

2 q > m + n + 2 : we know that 2 q > m + n + p, so either 2 q > m + n + 2 or 
2 q = to + n + 1 , in the latter since m + n + 1 < m + n+ p +1 < 2 g, we have 
p = 0 and this is covered by Case C. (3) “to + n + p + 1 > 0” becomes n > 0, 
which is evidently true by the definition of -D m ,n,p;q- The upper bound of the 
summation domain becomes [ n ~ m j. Therefore we can apply Lemma 

10 to the sum above and obtain 


D ri 


— Cn-l. n - 


q-l,n—q-\-l • 


5.2 m = 0 

When to = 0, similarly, for any n £ II vfc , 

7-)7T _ /-) n 02k+llT 

U 0,n,p;q — - Ly O, 2 fc I 0 ;fe+l- LA o,ra- 2 fe-l,p; l j-fc-l ' 

Again, summing over k and n £ II vfe we have 


Do,n,p;q — D m ,n,p;q — ^°’ 2k ’°’ k + 1 ^0,n-2k-l,p;q-k-l 

fc> 0 ien« k=o Trenvfc 

L^J LVJ 

= y \ T^o,2/c,0;fc-f-lT^o,n— 2k— l,p;q— k— 1 = y ^ Ck 1 kCq—k—2,n— q— k-\-l 


k -0 


k —0 


where the last equality comes from Case A, the induction assumption and the 
fact that Ck = Ck,k • Once again, we want to check the conditions in order to 
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apply Lemma 10 The condition “n > 0” is obvious, “to > p” is equivalent to 
2q > n + 3. Since 2q > n + p, there are two possibilities apart from “to > p”: 

1. 2q = n + p + 1. This is covered by Case B. 

2. 2q = n+p + 2 and p = 0. This is covered by Case A. 

“to + n + p + 1 > 0” is again equivalent to n > 0, which is evidently true. The 
upper bound of the summation domain “|^222±£J” j s Therefore we can 

apply Lemma 10 to the sum above and obtain: 

Do,n,p;q = Cq-l,n—q-\- 1- 


□ 


6 Proof of Theorem 7 

In this section we often abuse notations and do not differentiate between opera¬ 
tors and their kernels. Without loss of generality assume \v\ = 1 (otherwise one 
can scale (a, x , y,b)). We only consider the generating function of coefficient D , 
as the case for E can be dealt with similarly. We write 

m —1 n —1 

n,p)= ^2 Iaj (x)lA k ( y) II (xj-i - a - aA N ) (x k -i - Xj - /3A N ) 

1 <j<k<N ct— 0 /3=0 

p-1 

x Xk - 7An) 

7=0 

</>(to, n,p) = (x - a) m (y - x) n (b - y) p I a <x< v <b, 

where we recall A N = Then 0 < (f> N (m,n,p) < (j>(m,n,p ) < (b — a) m+n+p . 

We also write 


V 


m,n,p 


D 

J - y m,n,p;q / Q ^m+n+p+l- 2 q 

^ m\n\p\ 


Then 

JV-l 2JV-3 

w N = ^2 w s,n + 5 2 Vs ’ n 

s=l s=l 

where v s ,n are the degenerate terms, and 

w s ,n = 5Z 'Dr n ,n,p<t> N (.‘m,n,p). 
m-\-n-\-p—s— 1 

The reason for the range of the sum for s to be 1 < s < N — 1 is because 
w s ,n = 0 for s > N, as tj> N (to, n,p) = 0 for m + n + p > N — 1. 
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Proof of Theorem 7. We divide the proof into three parts: 

1- Em+n+p>jv «,p) N ^>°° 0 uniformly on [a, b) 2 . This shows 

the limit exists. 

2. w s,n converges weakly to W. 

3. The degenerate terms vanish uniformly: Y2 s =i 3 v s,n are arbitrarily small 
as N grows bigger. 


6.1 Part 1 

By the formula of D m ^ nyP - q a bound can be immediately obtained: 


0 <D m - q < ( Ln/2J ) <2". 

Similarly one can bound the trinomial coefficient ^ m +vn<p?' < 3 ™.+n+p. Combin¬ 
ing these two bounds we obtain 


|Pm,n,p| — 


m+n+p+1 


E 

9=0 


,D„ 


( 1 y l .m+n+p+l~2q J ^m,n,P 

m!n!p! 


< (m + n + p + 1) 


gm+n+p 

(m + n + p)! 


Therefore 


^ ' Pm,n,p0 (rU; Tlj p) 
m+n+p>iV 


< ^2 (m + n+p + 1 ) 

m+n+p> N 

V + 2 


— n\\ m + n +P 


(6(6- a)) 


(to + n + p)\ 


= E 


r>N 


(, + i) ll( ‘-. ,)l b o 


r! 


as AT —> 00 . 


6.2 Part 2 

We want to show that for any e > 0 and sufficiently large N, we have (/, X^ s <at-i( w; sW _ 
w s)g) < e||/|| 2 ||ff ||2 for testing functions f,g€ L 2 ([a, 6 )). 

Equivalently, we must show that 



f(x,y) 


E 

m-\-n-\-p<N —2 




(cj) N (m, n,p) - n,p))g(x, y)dxdy 


m,n,p 


<e||/|| 2 || 5 || 2 . 


We divide it into two further parts. 

1 . (x - a) m (y - x) n (b - y) p Ei<j<fc<jv J AixA fc (z, ?/) ~ (x - a) m (y - x) n (6 - 

2/) p Ia< x<y<b , 
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2. {x-a) m {y-x) n {b-y) p Y,\<j < k<N l A j xA k { x ^y) « Ei<j<fc<iv tN ti, k; m, n,p\ 
I Aj {x)I Ak {y); 

where 

m— 1 n— 1 p—1 

T N (j,k-,m,n,p) := (ajj_i - a - aA^) - Xj - /3A N ) JJ(6- x k - jA n ). 

a—0 (3=0 7=0 


6.2.1 Part 2.1 

We want to show that 

I" f f(x,y) Y V ^p E (x-a) m (y-x) n (b-yr 

Ja Ja m+n+p<N-2 l<j<k<N 


(a;)lA fc (y) - ^a<x<y<b)g(x, y)dxdy 


< e. 


Denote the left hand side by Bn, then by (1) and that (x — a) m (y — x) n (b—y) p < 
(■b — a) m+n+p for a < x < y < b, we have that 


B n < ~Y (m + n + p+l ) 

m+n+p<iV—2 


< (r + l) 


r<iV—1 



f(x)g(y)dxdy 

f(x)g{y)dxdy 


The term in the modulus can be bounded by repeated use of Cauchy-Schwartz 
inequality: 


N 

E 

1=1 


JV 


f(x)g(y)dxdy <Y \f( x )\ dx / lsG/)My 

Xj — i<x<y<Xj — i 

< E [ \f( x )\V x j~ x \\ghdx < Y Wfhhh = ^ Wfhhh- 


j =i 

Therefore 

iV-l 


i=i 


27V 2 


27V 


Siv < Y ( r + l)( r 2 2 ) m r[ a))r{b 21 f Wfhhh < CN-^WfUgh 

r =0 ' ' 


for some constant (7, where the second bound comes from the fact that X^ r >o( r + 

l)( r f)M^<oo. 
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6.2.2 Part 2.2 


We establish the following uniform convergence, from which weak convergence 
will follow: 

Y Vrn^p Y ^(aOlA*(y) (r N (j,k;m,n,p) - (x - a) m (y-x) n (b- y) p ) -> 0 

m-\-n-\-p<N — 2 l<j<k<N 

When a < x < y < b, t n is non-negative, and for m + n + p < N we use a 
telescoping series: 

0 - a) m (y - x) n (b - y) p - r N (j, k ; m, n,p) 

= (a; — Xj-i)... 

+ (xj -1 — a)(a; — Xj - 2 )... 

+ (xj-i - a)(a;j_2 - a)(x - £j_ 3 )... 

+ ... 

+ ( Xj-i - a)(xj -2 - a)(xj -3 - a)... ( x k+P -2 - y){b - y) 

+ (xj -1 - a)(a; i _2 - a)(xj_ 3 - a)... (6 - Xfc + p_ 2 )(a;fc+p-i - y) 

+ (xj_i - a)(a; i _2 - a)(xj_ 3 - a)... (6 - x fc+p _2)(6 - x fc+p _i) 

- (Xj-i - a)(a; i _2 - a)(xj_ 3 - a)... (6 - a; fc+p _ 2 )(6 - a: fc+p _i) 

( m n p 

E a +E^+E^ 

ct—1 /3=1 7=1 

< IV _1 (& — a) m+n+p ^(TO + n + p)(m + n +p + 1). 

Therefore 

E ^m,n,p E ^ ( X )^-A k {y) (t (j, k\m,n,p) 

m-\-n-\-p<N 1 <j<k<N 

-(x- a ) m (y - x) n {b - v y) < N- 1 Y ( k 2 2 ) + i) 2 (6(6 ~, a)) ' ^ CN ~ 1 - 

r<N ' ' 

6.3 Part 3 

The degenerate terms are the total orderings of path of length s with some 
repeated essential coordinates. If such a total ordering J has s + 1 — d non- 
repeated coordinates, then we call d the degree of degeneration, or we say that 
there are d degenerations in J. Each degeneration happens on a wedge part of 
a path, that is, any two essential coordinates Pi 1 j 1 = Pi 2 ,j 2 if and 0I1 ly if they 
correspond to parts of the same A fe part of a path for some k. On the other 
hand, degenerations happen in pairs. That is, for an array p , there do not exist 
three essential coordinates equal to each other, which would violate the partial 
ordering. Therefore, given a path of length s — 1, the number of degenerate total 
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orderings with d degenerations is bounded by 
paths of length s — 1 is the Fibonacci number 


8 -i 2 

J . Moreover, the number of 


$ s+2 _(_$)- 

v/5 


where $ = 


V5+1 


Since a path of length s — 1 can have at most wedges, there are at most 
^- degenerations. 

Therefore for each N, we have 


2JV-3 

Vn = ^ v s ,n 

s—1 


where V s< n is the counterpart of W s ,n that collects all degenerate cases of paths 
of lengths s — 1. By applying the calculation in the proof of Lemma 3 ( we can 
see the degenerate (of degree d) version of h™ is 


£{x,y) = 


(.b - a) d 
m\n\p\N d 


<t> {m,n,p) 


where m+n+p+d=s— 1. Thus the sum is bounded uniformly by 


(b-a) 8 - 1 

m\n\p\N d 


Therefore the total sum of degenerate terms is bounded as follows: 


|Vjv| < 


N-l V- 


EE E c ; *'« 


, (f> - a) 


8 -1 


s=l d— 1 m-\-n-\-p=s— 1—d 
N-l Er 1 

^ EE E cr 

s =1 d—1 m-\-n-\-p=s—l—d 


m\n\p\N d 


E(-d 


<?,,s- 2 9 


9=0 




( S -l)! 


sAT 


— d 


N-2 

— 2 — AT—1 

EE C 

d—l s=2d+2 


(6$(6 — a)) s 1 s^s —d + 2 
(3-1)! 


N — 2 , r , 

2 AT—1 


< E E < C( JV-. + £ «-) < cw-. 

d=l s=l ^ d=2 


□ 
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7 The unitarity of W 

Let / and g be the generating functions of D and E: 

m+n+p+1 

/:= E E D m ^ p . q [rn,n,p](-D) q u m+n+p+1 - q 

m,n,p>o g=0 

m+n+p+1 

g ■= E E - E m,n,p; 9 [m,«,p] t (-i /) <!, ^ m+ " +P+1_9 

m,n,p>o q —0 


E/-> 

S>1 

E&- 


And the kernel of W — I is 

f(x, y) < b a {x, y) + g{x, y ) (a:, y). (17) 

One can write down the equation that / and g have to satisfy for W to be 
unitary. 

Proposition 13. For W to be unitary, it suffices to show that for any a < x < 
V<b, 

_ rx _ ry _ rb _ 

f{x,y) + g(y,x)+ g(x,z)g(y,z)dz+ f(x,z)g{y,z)dz+ f(x,z)f(y,z) = 0. 
J a J x J y 

( 18 ) 

Proof. For W to be unitary it is necessary and sufficient to show it is both a 
coisometry and an isometry 

W*W = WW* = I 


Plugging in (17) and using the formulas for kernels of products and adjoints of 
integral operators we obtain equation (18) and three “other” equations: 

_ ry _ rx _ 

f(y,x) + g(x,y)+ g(x,z)g(y,z)dz+ g(x,z)f(y,z)dz 

J a J y 

+ [ f{x,z)f(y,z)dz = 0, a < y < x < b (19) 

J X 

_ r x _ ry _ 

9 {y,x) + f{x,y)+ f(z,y)f(z,x)dz+ f{z,y)g(z,x)dz 

J a J x 

+ / g(z, y)g(z, x)dz = 0, a < x < y < b 

J V 


_ ry _ r x _ 

f(y,x) + g(x,y)+ f{z,y)f(z,x)dz+ g{z,y)f{z,x)dz 
J a J y 

+ / g{z, y)g(z, x)dz = 0, a < y < x < b 

J X 


( 20 ) 


( 21 ) 
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The equation (18) and (19) are equivalent: one can interchange x with y and 
take a conjugate in the former to obtain the latter. So are (20) and (21). By 
the symmetry of / and g from Corollary 5, (18) and (20) are equivalent, hence 
it suffices to verify (18) to show the unitarity of W. □ 

Proof of Theorem 8. We divide the proof into two parts: first we write down / 
and g in a more amenable form, then we proceed to proving the integral identity. 


7.1 The formulas for / and g 

As previous, without loss of generality suppose \v\ = 1. Recall that 

„_, \n+j n+i n 

B A x iV) = E L , = {-lYWvy/tj&y/xy). 


n> 0 


(n + j)\n\ 


SillCG -£'rn.n.o:o — U-m— 


m,n,p;q — U-m=p=q,n= 0? 


g{x, y) = v^2 (-l) m —— ^ ^ = vB 0 (y - a,b - x). (22) 

z J mm 


m> 0 


Let f e and f Q be the parts of the sum of / where m + n + p are even and odd 
respectively, i.e. 

m+n+p+1 

f e = E E { x ~A m {y- x ) n {b-yY{-v)^ m+n+p+l - q 

2 | ra+n+p q —0 

ra+n+p+1 

fo= E E (x-a) m {y-xr(b-yn- 

2 fra+n+p q= 0 


i/yi/W+n+p+l-q 


The function f e can be further divided into the ( _? m ) part and the rest. Let 
u := (x — a)v, v := (y — x)v, w := (b — y)v and z = — ^ = —v 2 , then 


'E E 

g>0 m-\-n-\-p=2q 

= *E E 

g>0 m+n+p=2q 

= *E E 


u m v n w p z q ( n 


m\n\p\ \q — m 


w p z q f 2q — m — p 


m\(2q — m — p)\p\ \ q — m 

imv 2q-m-p w p z q / 2q _ m _ p 


q >0 0<m,p<g 
Q 


m\(2q — m — p)\p\ \ q — m 
{u/v) m ( w/v) m 


= vYv 2q z q E y _E 

Y'o ^ o m\(q-m)l^ o m\(q-m)\ 

(u + v) q (v+ w) q z q 
= -—-= vB 0 [y- a,b- x). 


q> 0 


q\q\ 


( 23 ) 
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The rest of f e is slightly more complicated. Let k = m+ " +p . We observe that 



Therefore the rest of f e , i.e. the sum excluding the terms corresponding to 

L-J * 


E 


u m v n w p 


E 


m\n\p\ 

k >0 m-\-n-\-p=2k q—k 


= v(z 
= u(z 
= u(z 
= v(z 
= v(z 
= v(z 
= v(z 


2k n 

dEEE 

k>0 n=k q=k 
2k n 

dEEE 

k ~>0 n=k q—k 
k k-\-n 

dEEE 

k >0 n=0 q—k 
k n 

dEEE 

/c>0 n=0 < 7=0 


n\ (it + w) 


2k—n yU 


qj (2 k — n)\ n\ 
z q (it + 


\2k—n yU 


q\ (2k — n)! (n — q)\ 


z q (u + io) 


k—r 


,,k+n 


q\ (k — n)\ (k + n — q)\ 

z k+q (it + -u>) fe - n t; fe+n 
(fc + g)! (k — n)\ (n — q)\ 


k k _ q ~k+q ( I , \k—q—n k+q+n 

dEEE* 1 ’ 


k >0 q=0 n =0 


(k + q)\ (k — q — n)\ n\ 


k 


DEE 

k >0 <7=0 


DEE 

< 7>0 fc>0 


(zi;) fc+9 (it + d + io) fc q 

(k + q)'- ( k-q)\ 

(zv) k+2q (u + v + w) k 
(k + 2q)\ k\ 


(24) 


We keep (24) to later merge it with a similar term in f 0 . 
For f 0 , let k = m +”+P~ 1 , We observe that 


9 


n z q 



Z 


k 



n 

k + 1 


y k +1 


+ ••• + 





Z 


n+1 


(*- 



Z 


fc+1 
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Following the same procedure which leads (24), the sum corresponding to the 
first term in the RHS is 


( zv) k+2q+1 (U + V + w) k 
y{ - z >2^2^ ( k + 2q + l)\ k\ 

q >0 fc>0 V H ’ 


(25) 


whereas the contribution from the second term is computed as follows: 


'£ £ 


u m v n w p (n 


z -—' ' m\n\p\ V k 

k >0 m+n+p=2fc+l x 

2fc+l 


-££ 


__ vfc+l 

\n\r>\ 

v n ('u + w) 2k+1 ~ n 


k>0 n=k 
k+1 


n\ (2/c + 1 — n)\ \k 

,=k x 


~k+l 


-££ 


,n-\-k 


(u + w) 


fc+1 —r 


(26) 


— J n\k\ (k + 1 — n)! 


~k+l 


k> 0n—0 

,.k ~/c+l 


n" z k+1 (v + u + w) k+1 
= v 2^ - 


k> 0 


k\ 


(k + 1)! 


= Bi(b - a,y- x). 


By summing up (23), (24), (25) (26) and plugging in u,v,w,z we obtain 
f(x, y ) = vB 0 (y - a, b - x) + B^b - a, y - x) - {is + v) ^ B g (y - x, b - a)v q . 

q> o 

(27) 

Note that E 9 >o B <i(y ~x,b~ a)u q = £ g > 0 M 2 V( b - a )iv ~ x )) (-y / £f 
is a generating function of the Bessel functions. 

7.2 Verifying the identity (18) 

We list a few useful properties of Bj (where we let B_\{x,y) := —Bi(y,x)): 

1. d x Bj(x,y) = —Bj_i(x,y),j > 0, 

2. B 0 {x,y ) = B 0 {y,x), 

3. Bj (0, y) = 6j 0 ,j > 0, 

4. d v Bj(x,y) = B j+ i(x,y). 

And an integral: 

fb 

/ B 0 (b - x,z - a)Bj{z -y,b- a)dz = - y~) B k {b -x,z- a)B k+j+1 (z -y,b-a) 
Jy k> n 


= f-Efe>o B k{b - x,b - a)B k+J+1 {b - y,b - a), j > 0 

{J2k>o B k( b - x ’ b - a) B k{b- y,b - a) - B 0 (b - x,y - a) j = -1. 


(28) 
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Substituting for / from (27) and g from (22) into (18) gives 

[ b 

Bi(b — ci,y — x)+ / Bq(z — a,b — x)Bq(z — a,b — y)dz 
J a 

f b 

+ / B\(b — a,z — x)Bi(b — a,z — y)dz 
■'y 

f b 

+ vBf){y — a, b — x) + vBq{x — a, b — y) + v / B\(b — a, z — x)Bq{z — a,b — y)dz 

J X 

f b 

+ v / B\(b — a ,2 — y)Bo(z — a,b — x)dz 

Jy 

f b 

— {v + V) / ^ Bi(b — a, z — x)B q (z — y,b— a)u q dz 

v q >o 

f b V 

— (v + i>) / B\(b — a, z — y)B q (z — x, b — a)v q dz 

J V _^n 

— a, b — y)B q {z — x, b — a)v q dz 

— a, 6 — x)B q (z — y,b — a)v q dz 


q> o 



9>0 


— (j/ + i/) B q (y — x,b — a)v q = 0. 

q> o 


(29) 


_^ ^ a;" 3 

Let Gj(x) = X)zc>o fc!(fc+i)! ’ f° r a > 0 the following two integrals hold: 

f Go(az)Go(/3z)dz = (a — /3)~ 1 (axGi(ax)Go(/3x) — / 3xGi(px)Go(ax )), 

Jo 

f Gi(w)Gi(w + a)dw = oT 1 (zGi(z)Gq(z + a) — (z + a)Gi(z + a)G 0 (z)) + Gi(a). 
Jo 

The first of these two integrals is the well-known Lommel’s integral, see e.g. 
Section 11 and 94 of Bow58 . The second integral written in terms of an 
indefinite integral of the Bessel functions is 

h/z( v ^ T? ) dw 

= /3~ 2 (wJi{w)J 0 (\/w 2 + /9 2 ) - V™ 2 + P 2 Ji ( V w 2 + P 2 ) Jo{w)^ ,P> 0. 

(30) 
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This is a special case of the so-called Sonine-Gegenbauer type integral (see e.g. 
page 415 of [Wat95 ). However, the authors have not found an explicit formula 
like the one on the right hand side of (30) in the literature. 

By these two integrals we have 


f b 

/ Bq{z — a 7 b — x)Bq(z — a,b — y)dz 

J a 


= ((y - x ) 1 ((& - y)Bi(b — a, b — y)B 0 (b - a,b - x) 
— (b — x)B\(b — a,b — x)Bo(b — a,b — y)) 


f b 

/ B\(b — a, z — x)Bi(b — a, z — y)dz 
J V 


= Bi{b~- a,y~ x) - (y - x) 1 ((&- y)B 1 (b- a,b- y)B 0 (b - a,b - x) 
— (b — x)Bi{b — a,b — x)Bo(b — a,b — y)). 


Therefore the first and the second lines of (29) vanishes. 

By the integral (28) (and interchanging x and y when necessary), the third 
and the fourth lines are reduced to their real part. 

Now if the real part of v is 0 then we are done. Otherwise by subtracting 
the first and the second lines and the imaginary part of the third and the fourth 
lines from (29), and dividing the remainder by the real part we simplify 
the integral identity into 


fb 

Bo(y — a,b — x) + Bq[x — a,b — y) + / B\{b — a, 2 — x)Bq(z — a,b — y)dz 

J X 
fb 

+ / Bi(b — a, z — y)B 0 (z — a, 6 — x)dz — 2 ^ B q (y — x, b — a)D q 

' v q> 0 

rb _^ 

— 2 / ^ Bi(b — a,z — x)B q (z — y,b — a)v q dz 

q> 0 

rb _^ 

— 2 / ^ Bx(b — a,z — y)B q (z — x, b — a)v q dz 

Xy q> 0 
fb 

— 2 v j ^ B 0 {z-a,b- y)B q (z — x,b — a)D q dz 

■' x q >o 
fb 

— 2v ^ B 0 (z — a, b — x)B q (z — y,b— a)v q dz 

■'y q> o 

rb _^ ^ 

+ 2 (u + D) / ^ B q (z — x, b — a)v q ^ B q (z — y,b — a)v q dz = 0 
Jy q >o 


33 


By the integral formulas (28), the above identity can be further simplified 
to 


E B k(b ~x,b- a)B k {b -y : b-a) + J2J2 B k (b - x, b - a)B k+q (b - y,b- a)i 

k> 0 g>0 k>0 

+ E E Bk ( b ~y^ b ~ a)B k+q (B -x,b- a)u q - ^ B q (y - x, b - a)v q 

q> 0 k> 0 q>0 

— / B\(b — a, z — x)B q (z — y,b — a)i' q dz 

' y g> 0 

f b \ 

— / Bi(b — a, z — y)B q (z — x 1 b — a)v q dz 
v q> o 

rb _^ _ 

+ (i/ + p) / B q (z — x, b — a)v q B q (z — y,b — a)u q dz = 0 


9>0 


q> 0 


It remains to verify the coefficient of v q for each g 6 Z in the LHS is 0, which 
can be done by repeated use of integration by parts and the properties of the 
Bj functions. 

When q = 0, the coefficient is 


J2 B j(b-x,b- 

- a)Bj(b — y,b — a) 

- B 0 (y - x,b-a ) 




3> 0 






+ [h-B^b-t 

J y 

i,z — y)B 0 (z — x,b 

- a) — Bi(b — a, z 

- x)B 0 (z — y,b — 

a) 


+ 5 2 B j+ 1 ^ Z ~ 

x 1 b — a)Bj(z — y,b 

~ a ) + ^2 Bj+i{z 

-y,b-a)Bj(z- 

■ x, b — 

a)dz 

j> 0 


j> o 




= '52 B j( b ~ x > 

b — a)Bj(b — y,b — 

a) ~ B o{y -x,b- 

a) 



j> 0 






+ E f B A* 

— x, b — a)Bj + 1 (z - 

-y,b-a)+ Bj(z - 

- y,b- a)B J+1 (z 

-z,b- 

- a)dz 

= J2 B j(b-x, 

b — a)Bj{b — y,b — 

a) - B 0 (y -x,b- 

a) 



j> o 






■ E B ^ z 

— x, b — a)Bj{z — y 

,b-a)\ b y = 0, 





j> o 


where in the fourth line we use integration by parts and the properties of Bj. 
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When q > 0, the coefficient is 


rb 

22 Bk{b - x,b — a)B k + q {b - y,b - a) + / —Si(6 — a,z — x)B q (z — y,b — a) 
k> o •'y 

+ 22 B + (z — x,b — a)Bj +q -i(z — y,b — a) + 22 Bj(z — x,b — a)Bj +q+ i(z — y,b — a)dz 


j> o 


i> o 

b 


= 22 Bk(b — x,b — a)Bk+ q (b — y,b — a) + / ^ Bj(z — x,b — a)B j+q - 1 (z -y,b-a) 

•’y j >o 


+ 22 Bj{z — x,b — a)Bj +q+ i(z — y,b — a)dz 


3 >~ 1 


r b 

222B k (b-x,b-a)B k+q (b-y,b-a)+ / 2^ Bj(z - x,b — a)B j+q - 1 (z — y,b- a) 
k> 0 ■’y j> 0 


+ 22 Bj-i{z — x, b — a)Bj +q (z — y,b — a)dz 


3 > 0 


= 22 Bk ( b — x ' b — a )Bk+ q {b — y,b - a) — 22 B k {z — x, b — a)B k + q (z — y,b — a)\ b y = 0. 


When q < 0, the coefficient is 

22 B k{b — y,b- a)B k+q (B - x,b - a) - B q (y - x,b- a) 

k> 0 

-Si(6 - a,z- y)B q (z - x,b - a) + 22 B k + q -i{z - x,b - a)B k {z -y,b- a) 

k> 0 

+ 22 Bk+ q +i(z — x,b — a)B k (z — y,b — a)dz 
k> o 

= 22 Bk ( b a)B k+q (B — x,b — a) — B q (y -x,b-a) 

k> 0 

+ 22 Bk+q~i(z — x,b — a)B k {z — y,b — a) + B k+q (z — x,b — a)B k ~i(z — y,b — a)dz 
k> o^y 

= 22 Bk ( b ~V’ b ~ a )B k+q {B — x,b — a) — B q (y -x,b-a) 

k> 0 

~22 Bk +<i( z - x,b - a)B k {z - y,b - a)\ b y =0. 

fc>0 



□ 
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